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VORTEX  FORMATION  IN  CHANNEL  BETWEEN  Ul'CKS  OF  FRU  ITON  THRUST  BEARINGS 

fArticle  by  M.  Ye.  Podol’skiv;  Moscow,  Mashinovedeniye,  Russian,  No  5, 

1970,  pp  92-99} 

The  load  capacity  of  friction  thrust  hearings  is  not  only  determined 
by  phenomena  that  occur  in  the  lubricant  fi ‘m,  but  also  depends  largely  on 
what  happens  in  the  bearing  as  a  whole.  Th.s  is  explained  basically  by 
three  ciicumstances.  First,  h  at  transfer  through  the  blocks  and  collar 
and  thermal  processes  in  the  film  that  *.«  transported  through  the  channel 
between  the  blocks  from  beneath  one  block  to  the  other  have  a  considerable 
influence  on  the  temperature  layer  of  the  film.  Since  the  viscosity  depends 
on  temperature,  the  noted  phenomena  have  a  direct  effect  on  the  magnitude  of 
the  hydrodynamic  reaction  of  the  lubric.  nt  wedge.  Second,  the  transfer  of 
heat  to  the  oil  flowing  in  the  body  of  the  bearing  is  accompanied  by  uneven 
heating  of  the  blocks,  which  causes  warping,  distortion  of  the  geometry  of 
the  lubricating  film  and  in  great  raensujc  leads  to  reduction  of  -its  load 
capacity.  Finally,  vacuum  zones  may  foi-ra  wjthir-  the  bearing  under  known 
conditions,  which,  as  shown  experimentally  !(•],  has  a  negative  effect  on 
working  capacity.  Analysis  of  the  block  cli  innei  is  obviously  very 
important  in  investigation  of  the  hydio.yn.-naics  and  beat  exchange  in  the 
body  of  a  hearing,  since  oil  flews  diicctly  from  it  into  the  oil  wedge. 
Experimental  studies  [2,  6]  and  also  on  determination  of  the  optimal 

filling  coefficient  [7]  show  that  pnn  in  the  block  channel  actually 

have  a  considerable  influence  on  hearing  operation  and,  in  particular,  on 
its  thermal  regime  (2).  Here,  so  far  -  on a  he  judged  according  to  [6}, 
the  main  feature  of  oil  motion  in  the  «. ‘tanne  1  is  the  formation  of  longitudinal 
vortexes.  The  latter  arc  created  as  le^nl*  of  rotation  of  the  volumes  of 
oil  between  the  blocks  under  the  infhw  n-e  of  forces  of  friction  from  the 
bearing  collar. 

An  attempt  is  made  below  to  earr.  out  an  approximate  analysis  of  the 
hydrodynamic  phenomena  that  e.cci.r  in  f  t-  block  channel,  which,  in  turn,  can 
serve  as  the  basis  for  refining  thermal  calculations  and  explaining 
conditions  that  cn  arc  the  absence  of  nenum  zones. 


The  problem  is  laid  out  schematically  as  follows.  The  channel  is 
replaced  by  a  recess  of  rectangular  croas  section  (Figure  1),  the  open 
siue  of  which  is  bounded  by  a  moving  plate.  The  distance  between  the 
plate  and  the  walls  of  the  groove  is  h  and  the  depth  of  the  groove  is  b. 
It  is  also  assumed  that  there  is  no  motion  of  the  lubricant  along  axis  z. 
Analysis  of  twisting  of  the  flow  with  consideration  of  longitudinal  flow 
along  the  channel  is  described  in  the  second  part  of  this  work. 


The  plane  problem  of  the  motion  in  the  slot  was  solved  [4,  11]  for 
h  =  0  and  IS]  for  h  t  0.  In  these  works  the  complete  system  of  Jfevier- 
Stokcs  equations  is  examined  and  they  are  numerically  integrated  with  the 
aid  of  computers.  Due  to  the  limited  memory  volume  of  the  computer,  however, 
the  solution  was  found  only  for  a  comparatively  small  number  of  Reynolds 
numbers  Re  =  Ub/v.  Meanwhile  the  problem  can  be  simplified  for  high 
Reynolds  numbers  (for  the  block  channel  Re  =  10**-105)  by  dividing  the  flow 
into  two  zones:  1)  boundary  layer  and  2)  nonviscous  core,  in  which  the 
forces  of  friction  can  be  ignored. 

In  this  statement  for  the  case  h  =  0,  as  it  applies  to  the  diagram 
in  Figure  2  (wall  1  is  fixed  and  wail  2  rotates  around  axis  0),  the  problem 
was  examined  in  [10,  12-14].  Here,  in  view  of  the  symmetry  of  flow  it  was 
assumed  that  the  pressure  gradient  along  the  boundary'  is  equal  tc  zero  and 
the  velocity  on  the  periphery  Ue  =  const. 


Figure  3. 
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On  !he  basis  of  analysis  of  the  equations  of  laminar  asymptotic 
boundary  1  -yer.  written  in  the  Prandtl-Miscs  form,  the  authors  of 
1 10,  12,  1  lj  showed  that  the  mean  square  velocity  in  the  boundary  layer 
along  each  His-i  is  the  same  and  equal  to 


VK  —  coast 


(1) 


Hence  for  velocity  U°  is  found  the  expression 


V'-Vcy s.  s  =  l/L, 


(2) 


where  Uq  is  the  velocity  of  wall  2,  It  was  also  shown  [10]  that  in  the 
general  case  of  plane  flow  with  closed  streamlines  the  vortex  rot  v  in  the 
nonviscous  core  has  a  constant  value. 

The  application  of  this  approach  to  cavities  of  rectangular  fora 
poses  known  difficulties  in  connection  with  the  fact  that  the  pressure 
through  the  length  of  the  boundary  layer  is  constant  and,  moreover,  vortex 
formation  is  possible  in  angles  A  and  B  of  the  grooves.  Angular  vortexes 
were  actually  discovered  both  as  a  result  of  numerical  solution  of  the 
complete  Navier-Stokes  equation  system  [4,  5,  11]  and  experimentally  [12]. 
Meanwhile  experimental  investigation,  carried  out  for  a  groove  of  square 
cross  section  for  Re  =  105  (the  pressure  here  was  laminar)  revealed  that 
the  pressure  gradient  along  the  perimeter  of  the  groove  is  small  and  the 
velocity  along  the  periphery  of  the  core  is  approximately  constant  and 
equal  to  U°  =  0.5U  [12].  This  result  can  be  explained  by  the  fact  that 

the  formation  of  angular  vortexes  leads  to  separation  of  the  streanlincs 
1'iom  the  walls,  with  the  result  that  their  form  approaches  circular  and 
the  flow  in  the  nonviscous  core  becomes  close  to  the  flow  according  to  the 
diagram  in  Figure  2.  As  regards  the  numerical  value  of  velocity  along  the 
periphery  of  the  vortex  Ue  =  0.5l»^,  it  can  be  found  from  (1)  for  L  =  4 1, 

i.e.,  schematic  representation  of  the  flow  in  the  form  of  a  flow  core  and 
unbroken  laminar  boundary  layer  yields  suitable  results  for  grooves  of 
square  cross  section.  Apparently  the  pressure  gradient  along  the  stream¬ 
lines  can  be  regarded  equal  to  zero  only  for  square  grooves  '.  In  grooves 
of  rectangular  cross  section  the  difference  of  the  streamlines  in  the  core 
from  circles  will  be  more  pronounced,  leading  to  variability  of  velocity 
and  pressure  along  the  periphery  of  the  core.  Meanwhile  in  the  practical 
sense  the  first  case  is  of  greatest  importance,  since  the  form  of  cross 
section  cf  the  block  channel  in  most  bearings  is  close  to  square.  From 
now  on,  therefore,  the  pressure  gradient  will  be  assumed  equal  to  zero. 

Examined  below,  in  the  specified  statement,  is  the  problem  of  vortex 
formation  in  the  channel  with  consideration  of  the  effect  of  tnc  lubricant, 
transferred  from  beneath  one  block  to  beneath  the  other  (h  j*  0) .  At  the 
end  of  the  article  arc  presented  certain  considerations  of  twisting  of  Jic 
flow  in  the  turbulent  flow  regime. 
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We  wi 1 1  examine  the  motion  of  the  lubricant  according  to  the  diagram 
in  Figure  3. 

H'e  will  separate  the  region  filled  with  oil  into  three  zones: 
nonviscous  core  III,  boundary  layer  II,  oil  stream  I,  forced  from  beneath 
the  preceding  block.  We  will  assume  that  there  is  a  separating  streamline, 
which  separates  zones  I  and  II.  From  the  physical  standpoint  this  assump¬ 
tion  is  quite  natural.  It  is  also  substantiated  by  the  results  of 
numerical  solution  [5].  Since  the  thickness  of  the  layer  in  zone  I  is 
small,  the  motion  of  the  lubricant  in  region  I,  as  m  region  II,  can  be 
described  by  the  boundary  layer  equations.  The  latter,  under  the  condi¬ 
tions  of  the  examined  problem,  are  conveniently  taken  in  the  Prandtl-Mises 
form  [3].  Then,  assuming  3p/3x  =  0,  we  obtain 

i  dn.'  d*u.' 

T  —  '  “  T“  • 

v  dx.  d\f." 

•  Here  x*  is  a  coordinate  on  the  streamline,  is  the  stream  function 
(u.  =  •1’/'*/3y* ,  v#  =  u*,  vA  are  the  longitudinal  and  lateral 

velocity  components,  respectively,  y„  is  the  lateral  coordinate. 

The  order  of  the  solution  of  the  problem  is  as  follows.  We  are 
given  an  expression  for  the  velocity  on  segment  abc,  separating  the  stream¬ 
lines.  I.et  this  expression  contain  n  unknown  parameters.  Then,  after 
determining  the  velocities  in  zones  I  and  II,  wc-  calculate  the  fcrce  of 
friction  on  line  abc  approaching  it  from  zones  I  and  II,  respectively. 
Having  now  broken  line  abc  down  into  n  segments  and  requiring  that  the 
forces  of  friction  on  each  of  them,  determined  dv  the  above-stated  method, 
art-  equal  to  each  other,  "we  find  n  equations  fer  determination  of  n 
unknown  parameters.  Further,  we  restrict  ourselves  to  the  simplest  case 
of  the  uniparanctric  problem.  Then  all  that  remains  is  to  calculate  the 
summary  force  of  friction  on  the  streamline  abc. 

To  find  the  approximate  solution  of  the  problem,  following  [12,  14], 
wc  will  linearize  equation  (1),  averaging  the  coefficient  in  front  of 
3?u*2/^*2. 

In  addition  we  will  proceed  to  the  dimensionless  variables  and 
introduce  the  Reynolds  number  Re,  assuming  (L  is  the  dimensional  length 
of  contour  abede) 


x 


2or. 

~T' 


Re  = 


UJL 


2av 


(3) 


Then  wc  obtain  approximately  for  zones  I  and  II 


da*  .  dV 
dz~*1  dfz' 


'f*  j  _  .  _  Utm 

UtL '  **  “  4 zrR'  Xl  _  4a*Re  ’ 


CD 


da* 

dx 


where  u,  and  are  the  average  valuer  of  the  dimensionless  velocity  in 
lm  2m 

zones  I  and  II. 

Satisfying  the  condition  of  periodicity  in  terms  of  x  in  zone  II, 
we  seek  u2  in  the  form  of  the  rea!  part  of  the  series 

or  -  ncut  (■•{»)  exp  (iA*x). 


Substituting  thi:.  into  the  -.econo  of  equations  (4),  we  obtain  for 
a^'i}/)  the  equation  system 


a  ak 


Viev'iig  *he  boundary  lever  as  a-.ymptotic  [10,  12-14],  we  find  that 
1  .  .*  vt  *2cit  ..  1  therefore  functions  a^vo  should  be  bounded  at  ip  00 . 

"olving  ti  -  .tions  for  a,  ;>i  wjrh  consideration  of  this  circumstanc-- , 
we  find  ? 


or  =  a.»  -f  2 sem. m.  ^7u«  \p[—  \k(l  -f  i)i\z  -rifcr],  rp=- 


ihc  coefficients  a,  „  am*  a,,  are  determined  as  a  function  of  the 
k0  C ') 

form  of  the  function  u(x)  on  the  separating  streamline  abc.  Ke  will 
consider  that  tp  =  0  on  line  abc  an  I  assume 


u'  —  d>cos- 


(*j  0,  —ns  <  z  -<  ns). 


Then,  .  ssuming  that  velocity  u  =  0  on  idea  and  expanding  u?(x,  <{>  = 
0}  into  a  Pourier  series  on  the  segment  -t>  <  x  <  r,  we  obtain 


a„  —  ait  -= - r  -  - 


<l>„  cos  kns 

ir  V  ~-(l:2*)s'  5 


om=— a». 

n 


The  summary  force  of  friction  on  the  separating  streamline  will  he 
fwith  an  accuracy  to  a  cons rani  cofactor) 

J  f  du1  | 

T  f -  dx. 

-  id'H-  1  «»• 


(8) 


Substituting  here  u2  fiom  (5)  and  a,n  from  (7),  we  find 


T  = 


Xt 


kO 

m.  /<«>=-4 Y-J^t 

Y 2~}rk(l-4ktsi) 


/<°,25)  =  -l{£ 


(-1)’ 


f9) 


JLl 


V2n  +  1[4 -(2r  -f  1)*]  8f2  * 


ri}- «  0,5. 


We  will  now  examine  zone  I  (the  first  of  equations  (4)).  Equation 
( i  j  must  be  solved  for  the  conditions 


ir=:— ^(D,,  -os— -  *  for  «t==0,  u*==i  for  $  =  —  c.  (10) 

2  2s 

Here  (-c)  is  the  stream  function  on  the  moving  plate.  If  we  assume 
that  the  oil  flow  rate  through  unit  width  of  zone  I  is  0.5U_h,  then  (sec 
(1))  c  =  O.Sh/L.  u 

We  will  seek  the  solution  in  the  form 

a1  =  *  +  a,ft>)exp  J-  a-,(ij>)exp  (11) 


where  in  view  of  (10)  the  functions  (ij>)  and  a  j(i|»)  should  satisfy  the 
conditions 


(0)  =  (0)  =  —  <Vm,  o«( — 6,)  =  a,(— 6,)  =  0,  6, 

4 


2xt*  '  W 


Substituting  (11)  into  (4)  we  find  for  and  a  ^  differential 

equations  of  tne  second  order.  Solving  these  equations  for  conditions  (12) 
and  then  substituting  otj  and  a  ^  into  (11)  we  find 


u'~ - —  r  -rp  Bern.  *». 


I  sbl(^i  -0} 

sfa  16(l~z)j 


l 


2xj  Vs 


(13) 


Hence  for  the  force  of  friction  on  line  abc  on  approach  from  zone  I 
we  obtain  with  the  aid  of  (8) 


z(M  = 


7= - ^i-5.X(6.)3V}, 

(sh26,  -f  sin26,)(ch26,-f  i  s26,) 

sh1 26,  -f  sin:  26,  ~  ^ 
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Further  we  will  examine  the  case  of  small  6,.  Calculations  show  that 


when  6^  <  1  the  function  x(<$j) 


fij"1 ,  so  that 


-(i  — <!>„). 


Equaling  (9)  and  (14)  we  obtain  the  equation 


T =-£«i=  ns) 

n$x2  s  L 1  u;m 

2  2 
The  value  U2m  ,  on  the  basis  of  (1),  is  assumed  to  be  equal  to  u^ 

on  the  outer  boundary  of  the  boundary  layer  in  zone  II  [12,  14).  We  obtain 
(see  (5)  and  (7)) 

(16) 

Assuming  now  (for  a  square  groove)  s  =  0.2S,  we  find  from  (15)  and 


1  -  <H-  =  <p*d)„  ■%  <p°  =  2,83  - 


We  will  assume  h  =  20  mkm  =  20*10-6  m,  =  50  m/sec,  v  =  O^'IO”1* 

m?/soc,  L  =  8*10-2  m  (to-  this  value  of  L  corresponds  a  block  channel  width 
of  approximate1}’  2  cm).  Then  R  =  3.2*10"'‘.  <J>°  =  0.13  and  equation  (15) 
yields  <I>  -  0.9\  The  dimensionless  velocity  on  the  periphery  of  the 

®  a  ^  A  ^ 

vortex,  in  view  of  (15)  is  U  =  /a__  =  /s<P  or  for  s  =  0.25  U  =  0.5/F*. 

’  '  00  m  m 

It.  follows  from  the  above  calculations  that  the  lubricating  film  decreases 

the  velocity  on  the  periphery  of  the  vortex,  but  this  reduction  is  small 

when  the  thickness  of  the  fill.',  is  small  (in  the  examined  case  we  have 

0  0 

instead  of  U  =  0.5  u'  =  0. 5/0.9  =  0.47).  This  finding  means  that  the 
flow  is  close  to  that  obtained  if  the  layer  is  absent  altogether,  i.e., 
if  h  =  0.  For  h  =  0,  however,  the  velocity  on  the  line  ip  =  0  -  vs  <  x  < 

<  ns  is  constant.  Therefore  in  the  examined  conditions  for  velocity  on' 

Rising  these  results  for  estimating  the  value  6^  from  (12)  and  assuming 
u  =  /0.5(1  +'  <»>  )  =  1 ,  we  obtain 

m«  v  m  9 


6,  =  —  2.7  )«  a  0.1S. 
21 


i  '  '  '  1 

*  •  •  I  *  1  J  • 
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the  separating  streamline  it  is  advisable  to  assume  u  =  =  const  instead 

of  dependence  (6).  finding  for  this  case  the  solution  analogous  to  the 
preceding,  we  again  obtain  for  4>  equation  (15),  but  now  f(s)  is  deter¬ 
mined  according  to  the  equation  m 

..  .  2  n  siir&ns 

/($)  =  -=->  - — — , 

/ (0,23)  “  ~y  i  +  i)  {  1  -  -U  5  (-)  -  0,649, 
nV?'  V2  /  V  2V2 '  \  2  / 


where  C(x)  -  £  is  the  Riemann  function,  C(-|)  =  2.61  [9]. 

For  the  above-examined  case  we  have  =  0.17,  4>  =  0.85,  u^  = 

' _  Q  ® 

=  0.5/0.85  =  0.46,  i.e.  u  close  to  the  value  obtained  above. 

We  will  now  examine  the  problem  of  vortex  formation  in  a  groove  in 
the  assumption  that  the  motion  is  turbulent.  Here  we  will  assume  h  =  0. 
In  this  case  the  region  occupied  by  oil  is  broken  down  into  two  parts: 
the  flow  core  and  turbulent  boundary  layer. 


The  boundary  layer  equations  will  be 

da.  da.  1  dx.  du»  dv» 

a.~  +  y.  — =- - r— ,  - - Ft— 

dx.  dy.  p  Oy.  dx»  uy» 


The  force  of  friction  t*  is  given  by  the  dependences  [1]  . 

=  n—  for  0^y.  sg6„  t  =  pk:y.2 (  for  (19) 

dy.  \dy./ 


where 


—  I  =t—'i 
d,J- 


Here  Sq  is  the  thickness  of  the  laminar  film,  p  is  the  density  of 
the  lubricant,  V  and  kj  are  constants. 

Writing  (18)  in  Prandtl-Mises  variables,  we  obtain 

da.  dx. 

Pt— =  — — :  (2f 

dx.  dip.  v 
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Integrating  (20)  in  terms  of  x*  within  the  limits  from  zero  to  L 
and  considering  that  on  the  outer  boundary  of  the  boundary  layer  friction 
is  equal  to  zero,  we  will  have 


4* 

j  t#  |  ».«o  dx%  =  0. 


Equation  (21)  represents  the  obvious  condition  of  equilibrium  of 
rlie  boundary  layer. 

Through  and  t ^  we  will  denote  the  average  forces  of  friction 

on  the  wall  on  two  segments, i respectively:  x.  e  (0,  l)  and  *.e{r,L  —  l) 

Then  we  have  from  (20)  [see  also  (2)) 

! 

=  —  s).  (22) 

•If,  following  Prandtl's  example  [8],  we  assume  that  the  force  of 
friction  is  constant  on  a  line  perpendicular  to  the  wall,  then,  by 
integrating  equations  (19)  for  the  segments  z.  e(0,/)  and  e  (/,/,—  /), 
we  obtain 

_ p£V(l-uy  pW 

t., - - - c„  =  - c,  ■  .  (23) 

where  Cj  and  c 2  are  the  coefficients  of  resistance,  determined  by  the 
formulas  (6}  and  are' the  boundary  layer  thicknesses) 

c,=  (— Jn— +— r!  6,  (i  =  1,2). 

'  k  S„  k  I  •  k  Vt.,  . 

Formulas  (23)  shew  that  the  ref .stance  of  friction  can  be  determined 
according  to  the  formulas  for  tubes  if  the  difference  of  velocities  on  the 
cuter  boundary  of  the  boundary  layer  and  on  the  wall  is  regarded  as  the 
characteristic  velocity.  Using  on  this  basis  the  interpolation  formulas 
for  the  coefficients  of  resistance  [1],  we  obtain 

„  =  0,225  C!  =  0io2,5  (24, 


We  will  now  consider  that  on  each  of  the  segments  the  thicknesses 
of  the  boundary  layers  are  approximately  constant.  Then,  from  (22),  (23) 
and  (24)  we  obtain 


(  l~ u°  M*  u'  _]/  * 

\  u®  6 ,/  1  —  u®  I  l-( 


KTD-  i  1 C-2  3  - 1 8  j  9  -  7 1 
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■  r  - 1 


The  ratio  of  thicknesses  6^/6^  can  he  found  by  assuming  that  the  oil 

flow  rates  in  films  1  and  2  are  equal  to  each  other.  Then,  assuming  for 

velocity  the  "one'Seventh"  wall  [2], 


u • 


we  have  =  7u°/(l  +  7u°)  and 

formula  (24)  yields1 

Ta(a‘) 

5  1+<P*(“#) 

,  r  7(1  — u#)  i’*  a8 


(26) 


‘  The  graph  of  the  function 
u°  -*  u°(s),  plotted  according  to 

formula  (26),  is  illustrated  in  Figure  4  (curve  2).  Also  shown  here  is 
curve  1 ,  plotted  according  to  equation  (2) ,  obtained  for  the  case  of  a 
laminar  boundary  layer.  Comparison  cf  curves  1  and  2  shows  that  the 
intensity  of  twisting  of  the  flow  core  in  the  turbulent  regime  is  less  as 
a  whole  than  in  the  laminar  regime.  For  s  =  0.25,  in  particular,  u°  =  0.36. 


The  twisting  of  the  flow  can  be  judged  indirectly  on  the  basis  of 
the  pressure  drop,  which  forms  through  the  cross  section  of  the  channel 
as  a  result  of  centrifugal  forces.  It  can  be  determined  according  to  the 
formula  Ap  =  O.SpVg2,  where  vQ  is  the  velocity  on  the  periphery  of  the 
twirted  core. 


The  values  (0.2-0.4)u^  are  given  in  [6]  as  the  approximate  values 

for  Vq.  The  latter  in  the  laminar  regime  of  motion  are  less  than  those 

that  follow  from  the  equations  derived  above  (for  channels  of  square 
form  vft  *  0.5U0  in  laminar  flow  and  vQ  a  0.36U0  in  turbulent  flow), 

which  is  explained,  in  particular,  by  the  fact  that  the  twisting  of  the 
oil  entering  the  channel  develops  gradually.  Therefore  the*  value 
Ap  =  0.5p(u“U0)2  should  be  viewed  as  the  limiting  value,  established 

during  complete  twisting  of  the  flow.  Analysis  of  this  problem  will  be 
given  in  the  second  part  of  this  work. 
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